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Abstract 

Motivated by an application to LDPC (low density parity check) algebraic geometry 
codes described by Voloch and Zarzar, we describe a computational procedm'e for 
establishing an upper bound on the arithmetic or geometric Picard number of a smooth 
projective surface over a finite field, by computing the Frobenius action on p-adic 
cohomology to a small degree of p-adic accuracy. We have implemented this procedure 
in Magma; using this implementation, we exhibit several examples, such as smooth 
quartics over F2 and F3 with arithmetic Picard number 1, and a smooth quintic over F2 
with geometric Picard number 1. We also produce some examples of smooth quartics 
with geometric Picard number 2, which by a construction of van Luijk also have trivial 
geometric automorphism group. 

Introduction 

Much recent work has gone into the computational problem of computing the zeta function 
of a "random" curve over a finite field, in large part because the question of determining the 
order of the Jacobian group (class group) of such a curve stands in the way of using said 
group for public key cryptography. The history of this problem is not our present concern, 
and anyway it has been documented elsewhere; see for example [27] for an overview. 

By contrast, relatively little work has gone into the analogous computational problem 
over higher dimensional varieties. Ongoing work of Bas Edixhoven and his collaborators, 
to give an efficient algorithm for computing the n-th Fourier coefficient of a fixed modular 
form when n is a large integer of known factorization, amounts to computing factors of 
zeta functions of higher- dimensional varieties over large prime fields using methods of £- 
adic cohomology. Over fields of small characteristic, one also may use techniques of p-adic 
cohomology, which when applicable tend to yield more efficient algorithms. However, while a 
number of reasonable-looking algorithms for various higher-dimensional varieties have been 
described theoretically, e.g., by Gerkmann [H], Lauder [30l[3T], and Lauder and Wan [33] , 
until recently nothing had been attempted in practice. (For some very recent developments 
on this front, see Section [3^ ) 
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The purpose of this paper is to begin repairs on this gap in knowledge, by on one hand 
illustrating how even limited information about the action of Frobenius in p-adic cohomology 
can be used to address questions of some possibly practical import, and on the other hand to 
outline an algorithm which has been demonstrated in practice to be able to obtain this limited 
information. The potential import stems from the fact that one can use information about 
Frobenius, specifically bounds on Picard numbers obtained from performing linear algebra on 
a low-precision Frobenius matrix, to control the minimum distance of an algebraic geometry 
(Goppa) code derived from a surface. As observed by Voloch and Zarzar, such codes have 
the LDPC (low density parity check) property and so may be of special interest. 

Besides this introduction, the paper is structured in four main sections. The first section 
is general, describing in detail what a Picard number is and how to use an approximately 
computed Frobenius matrix to bound it. The second part gathers some facts about algebraic 
de Rham cohomology and p-adic cohomology. The third part sketches a particular algorithm 
for producing an approximate Frobenius matrix on the cohomology of a smooth hypersurface, 
using p-adic cohomology and a description of the cohomology of a smooth hypersurface due 
to Griffiths [17J; we also mention some related proposals. The fourth part describes an 
implementation of our algorithm in MAGMA and tallies a few experimental results. 



1 Picard numbers and Frobenius matrices 
1.1 Picard groups 

Definition 1.1.1. Let X be a variety over a field k. The Picard group Pic(X) is the group 
of isomorphism classes of line bundles (or invertible sheaves) on X. Note that for X smooth, 
isomorphism classes of line bundles are in natural bijection with rational equivalence classes 
of (Weil or Cartier) divisors on X. 

Lemma 1.1.2. Let X be a smooth proper irreducible variety over a field k, let k^^^ denote 
the separable closure of k, and put G = GaA{k^'^^ / k) . Then the natural map 

Pic(X) ^Pic(X XfcF^P)^ 

(in which the superscripted G means take G-invariants) is always infective; moreover, if the 
Brauer group Br(fc) is trivial (e.g., if k is finite), then the map is surjective. 

Proof. For any smooth X, the Hochschild-Serre spectral sequence in etale cohomology gives 
rise to an exact sequence 

^ H\G, k'^P[X]*) ^ Pic(X) Pic(X Xfc k'^'Pf ^ H\G, fc^"P[X]*) 

where /i;''''P[X]* = H^{X k^'^^^Qrn) (see, e.g., [HI Lemme 6.3]). For X proper irreducible, 
k^^^[X]* = (/t'^^P)*, and H\G, (fc'^^P)*) = by Hilbert's Theorem 90, while H^{G, (/t'^^P)*) = 
Br(/c). This yields the desired results. □ 
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Definition 1.1.3. For X smooth proper irreducible over a field k, we say that a divisor is 
algebraically equivalent to zero if it has the form 

(Dn(XxM))-(Dn(Xx{g})) 

for some connected (but not necessarily smooth or irreducible) curve C, some pair of closed 
points p,q on C of the same degree, and some divisor D on X x C containing no fibres of 
the projection X x C C. The set of divisors algebraically equivalent to zero is a subgroup 
closed under rational equivalence; let Pic°(X) denote the image of this subgroup in Pic(X). 

Remark 1.1.4. For k algebraically closed, the elements of Pic°(X) can be thought of in a 
natural way as the /c-valued points on a certain variety over k, the Picard variety associated 
to X] there is also a scheme-theoretic version of this fact that works over more general bases. 
We will not use this interpretation here. 

Lemma 1.1.5. Let X he a smooth irreducible complete intersection in P^, for k a perfect field 
and r a positive integer. Then Pic°(X) = 0, Pic(X) is torsion-free, and 0{1) is indivisible 
zn Fic{X). 

Proof. For k algebraically closed, this is [lOl Theoreme 1.8]; the general case follows from 
the algebraically closed case plus Lemma 11.1.21 Note that if dim(X) > 3, one in fact has 
Fic{X) = Z-0{1). □ 

1.2 Neron- Sever i groups 

Definition 1.2.1. For X smooth proper irreducible over a field k, the quotient NS(X) = 
Pic(X)/ Pic°(X) is called the Neron-Severi group of X. 

Remark 1.2.2. For X projective, one may define the degree of a divisor with respect to 
any fixed ample divisor. The resulting map induces a homomorphism deg : NS(X) —>■ Z 
sending any ample divisor to a positive integer; in particular, NS(X) is nontrivial and any 
ample divisor represents a nonzero class in NS(X). 

Remark 1.2.3. For A; = C, there is a natural map 

NS(X) ^ i/2(X^°, Z) n H^'\X), (1.2.3.1) 

where X^^ denotes the analytic space associated to X and H^'^{X) = H^{X, fix) (which by 
GAGA may be computed either algebraically or analytically); the Lefschetz (l,l)-theorem 
[T8| §1.2] asserts that the map fll.2.3.ip is a bijection. 

Definition 1.2.4. Let X be a smooth proper irreducible variety over a field k. From the 
Lefschetz theorem, it follows that NS(X) is finitely generated in case k = C In fact, NS(X) 
is always finitely generated; this was first shown by Neron [39l Theoreme 2, p. 145] (see also 
[2Tj . [Ml §1]). The rank of NS(X) as a finitely generated abelian group is called the Picard 
number (or arithmetic Picard number) of X. The rank of NS(X x^k), for k the algebraic 
closure of /c, is called the geometric Picard number of X. 
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Definition 1.2.5. Let X be a smooth projective surface over a field k. Then intersection 
theory Chapter V] gives rise to a symmetric pairing on divisors of X, called the intersec- 
tion pairing; this pairing respects algebraic equivalence, so it descends to NS(X). We say a 
divisor D is numerically equivalent to zero ii C -D = for every projective curve C contained 
in X] this turns out to happen if and only if some multiple of D is algebraically equivalent to 
zero p6l Theorem 4]. That is, the quotient of NS(X) by the classes of divisors numerically 
equivalent to zero is precisely NS(X)/ NS(X)tors- This group is never zero because an ample 
divisor D satisfies C ■ D > for any C, and so is not numerically equivalent to zero. That 
is, the Picard number of a smooth projective surface is always positive. 

Remark 1.2.6. Although we will only attempt to bound Picard numbers over finite fields, 
doing so also has consequences over number fields. That is because if p is a prime ideal in 
the ring of integers Ok of a number field K, k = Ok/P, and X is a smooth projective surface 
over the localization of Ok at p, then the torsion-free quotient of NS(Xi^) injects into the 
torsion- free quotient of NS(Xfc), compatibly with the intersection pairing [49], §6]. We can 
thus control the size of NS(X/^) by controlling NS(Xfc); in some cases, one can gain further 
control by reducing modulo more than one prime of good reduction [50] . 

1.3 Picard numbers and codes 

We now recall briefly what Picard numbers have to do with error- correcting codes; the link 
lies in a higher-dimensional version of Goppa's construction [THj of algebraic geometry codes 
from curves over finite fields. 

Definition 1.3.1. Let X be a smooth projective irreducible variety over a finite field Fg. 
Let H be an ample divisor on X, let m be a positive integer such that the divisor mH is very 
ample, and let C{mH) = T{X, 0{mH)) be the Riemann-Roch space of mH\ we may identify 
elements of C{mH) with rational functions / G Fg(X) such that the divisor div(/) +mH is 
effective. Let S be the set of F^-rational points of X \ if. Define the code C{X, mH) to be 
the subspace of of functions induced by elements of C{mH), viewed as a linear code over 
F 

In Goppa's original construction, X is a curve, the rate of the code (the ratio between 
the dimensions of the code and of its ambient space) is determined by Riemann-Roch, and 
a good bound on the minimum distance (the smallest number of nonzero elements in a 
nonzero codeword) comes from the fact that a rational function cannot have more zeros than 
poles. In higher dimensions, one can still get rate information out of Riemann-Roch, but 
bounding the minimum distance is trickier. For surfaces, this question has been investigated 
by Voloch and Zarzar, with the idea of using the subcodes induced by curves on a surface 
to give an asynchronous decoding algorithm in the style of Luby-Mitzenmacher [35]. Voloch 
and Zarzar observe that a low Picard number gives rise to a good bound on the minimum 
distance; we limit ourselves here to mentioning two sample assertions, and defer to [51] for 
more information. 
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Lemma 1.3.2 (Voloch). Let X be a smooth projective surface over a field k, and suppose 
NS(X)/ NStors(-^) is generated by the ample divisor H . Then for any positive integer m, the 
zero divisor of a nonzero element of C{mH) has at most m irreducible components. 

Proof. See |5T1 Lemma 2.2]. □ 

Lemma 1.3.3 (Zarzar). Let X be a smooth surface of degree d in P'^ over a perfect field k, 
and suppose that the Picard number of X is equal to 1. IfYis an irreducible surface in 
of degree m < d, then X P[Y is also irreducible. 

Proof. First note that Pic(X) and NS(X) coincide and are torsion-free by Lemma ILL5I 
Then invoke [53l Lemma 2.1]. □ 



1.4 Zeta functions and Picard numbers 

Definition 1.4.1. Let X be a smooth proper variety over a finite field ¥q. The zeta function 
of X is the power series 

(oo rpn\ 
^#X(F,n)— j . 

The Weil conjectures, proved by Dwork, Grothendieck, Deligne, et al. (see [22l Appendix C] 
for a fuller statement), assert that there exists a product decomposition 

2dim{X) 
i=0 

where Pi{T) G Z[T] and Pi{0) = 1, such that the roots of Pi{T) in C all have absolute value 
g-V2_ Moreover, the multiset of roots of Pi is invariant under the transformation r i— > q~^/r. 

The connection between zeta functions and Picard numbers was first articulated by Tate 
[H], who showed that 

rankNS(X) < ordT=i/g P2(T). (1.4.1.1) 

(Actually Tate's argument gives a bit more information than this; see Remark 11.5.21 below.) 
Tate conjectured further (by analogy with the conjecture of Birch and Swinnerton-Dyer) 
that equality always holds in fll.l.l.ip : Tate himself showed this for abelian varieties, and 
it is also known in some other cases, e.g., for ordinary K3 surfaces [52]. Tate's conjecture 
in general is wide open; however, since our purpose here is merely to give upper bounds for 
Picard numbers, the unconditional inequality (11.4. 1.11) will suffice. (See |15] for more context 
on Tate's conjecture.) 

Note that (11. 4. 1.11) also gives a bound on the geometric Picard number: 

rankNS(X Xp^ F^) < ^ ordT=c/, P2(r), (1.4.1.2) 

C 
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where ( runs over all roots of unity. Since P2{T) has integer coefficients, we can rewrite 
f ll.4.1.21) as 

rankNS(X Xp, F^) < ^ 0(n) ordT=c„/, ^(T), 

n:</)(n)<dog(P2) 

where (n denotes any one primitive n-th root of unity. For computational purposes, we need 
an explicit bound on n; here is an easy such bound. 

Lemma 1.4.2. For any positive integer n, we have 

n 



n) > 



[log2(n)J +1' 

Proof. We have 

n \ p 

p ^ 

the product running over the distinct prime divisors of n. There are at most Llog2('^)J such 
divisors, so 

^ > TT = - 

as desired. □ 

Here is a standard parity consideration that arises in the context of Tate's conjecture. 

Remark 1.4.3. Note that the right side of (11.4. 1.21) has the same parity as ordr=i/g P2(T) + 
OTdT=-i/q P2{T). This in turn has the same parity as deg(P2)5 since il/g are the only real 
roots consistent with the restriction that each root has absolute value 1/q. Under Tate's 
conjecture, equality would hold in fll.4.1.2p . and would thus imply that if X is a smooth 
proper variety for which deg(P2) is even, then the geometric Picard number of X is at least 
2. 



1.5 Weil cohomologies 

Definition 1.5.1. Fix a finite field Fg of characteristic p, and let be a field of characteristic 
zero. By a (weak) Weil cohomology over K, we will mean the following data. 

• One must specify a collection of contravariant functors from smooth proper varieties 
X over ¥q to finite dimensional i^- vector spaces H^{X) equipped with endomorphisms 
Fj, such that 

PiiT) = det(l - TFi, H\X)) (i = 0, . . . , 2 dim(X)). 

For m G Z, we write H^{X){m) to mean the vector space W{X) equipped with the 
endomorphism q~^Fi. 
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• The Lefschetz trace formula holds for Frobenius: for any positive integer n, 

2dim(X) 

#X(F,.)= Yl (-lrtrace(i^^i^^(X)). 

• One must specify functorial, F-equivariant maps (for d = dim(X)) 

tracex : H^\X){d) ^ K 

(where F acts as the identity on K) which should be isomorphisms when X is geomet- 
rically irreducible. 

• One must specify associative, functorial, F-equivariant cup product pairings U : H^{X) x 
W{X) W+\X) such that (for d = dim(X)) the pairings 

W{X) X H^'^-\X){d) ^ H^\X){d) K 
are perfect (Poincare duality). 

• One must specify an injective i^-linear homomorphism 

NS(X) ®j^K^H\X){lf=^ 

(the cycle class map). 

For a more precise definition of the phrase "Weil cohomology" (which actually includes more 
structure than this, including a Kiinneth decomposition, cycle class maps for higher Chow 
groups, and a full Lefschetz hyperplane theorem, plus additional compatibilities), see [29] . 

Remark 1.5.2. By virtue of the cycle class map, the existence of a Weil cohomology yields 
the inequality (11. 4. 1.11) . as the right side of (11.4.1.11) equals the dimension of the generalized 
eigenspace of H'^{X) with eigenvalue q. In practice, we will use the resulting slightly stronger 
version of (11.4.1.11) : 

rankNS(X) < corank(F2 - q,H^{X)) (1.5.2.1) 
and the corresponding version of (11.4.1.21) : 

rankNS(X Xp^F^) < ^ 0(n) corank(F2 - C„g, i^^(X)). (1.5.2.2) 

?i:0(n)<deg(P2) 

In theory, one expects that F2 acts semisimply on if^(X); this would follow from the full 
conjecture of Tate, which is somewhat stronger than we have described here (as it makes 
predictions about if^*(X) for all i). 

Remark 1.5.3. At the time |44] was written, the only Weil cohomologies that had been 
constructed were the £-adic etale cohomologies for each £ p, which take values in Q^. 
Subsequently, it was shown that Berthelot's rigid cohomology is also a Weil cohomology; it 
takes values in the p-adic field Qq. (Here and throughout, for brevity we write Zg for W(¥q) 
and Qq for FracZ^.) For the essential properties of rigid cohomology, see [21 H]; also see [23] 
for additional context on j9-adic cohomology theories. 
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1.6 Approximate Gaussian elimination 

Let us set some notation for this subsection. 

Notation 1.6.1. Let K he a. complete discretely valued field. Let Ok be the ring of integers 
of K, fix a uniformizer vr of Ok, and write v{x) for the valuation of x & K. 

We are going to describe an algorithm for producing an upper bound on the corank of a 
matrix A over K, given only the information of the entries of A modulo vr™ for some integer 
m. There is no harm in rescaling the matrix A (by multiplying by an appropriate power vr" 
of TT, then replacing m by m + n) to ensure that A has entries in Ok- 

Algorithm 1.6.2. Given a matrix A with entries in Ok/tt^^Ok which is the reduction of a 
matrix A over Ok, the following algorithm returns an upper bound on corank(A). 

1. Let r be the number of rows of A, and let s be the number of columns of A. If A has 
no nonzero entries (possibly because A is an empty matrix), return s and STOP. 

2. Choose a nonzero entry Aij of minimum valuation. 

3. For k = 1, . . . ,r in succession, skipping over k = i, choose c G Ok/i^"'Ok such that 
cAij = Akj, and subtract c times the i-th row of A from the k-th row of A. 

4. For i = 1, . . . , s in succession, skipping over i = j , choose c G Ok/tt"'Ok such that 
cAij = Ai£, and subtract c times the j-th column of A from the i-th column of A. 

5. Delete row i and column j from A, then go to step 1. 

Proof. We prove the claim by induction on the number of rows plus columns of A. If A is the 
zero matrix, the claim is evident. Otherwise, we may assume without loss of generality that 
i = j = 1. By lifting each of the row and column operations from A to A appropriately, we 
may also assume that An = for i = 2, . . . ,m and that Ay = for j = 2, . . . , n. It is now 
clear that the corank of A is equal to that of its lower right (m — 1) x (n — 1) sub matrix. □ 

One can also use Algorithm IL6. 21 to obtain information about the determinant of A. 

Proposition 1.6.3. Given annxn matrix A with entries in Ok/t^^Ok which is the reduction 
of a matrix A over Ok, suppose that Alaorithm \1.6.2\ returns the bound 0. For h = 1, . . . ,n, 
with notation as in the h-th iteration of the algorithm, put Oh = A^j G Ok/t^"^Ok and Ch = 
(—!)*+■' G Ok- Ghoose lifts ai, . . . , a„ of to A. Then 

n 

{m-v{ai)} + } v{ai). (1.6.3.1) 

i=l 
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Proof. Perform the "shadow" computation of the proof of Algorithm II. 6. 21 with the following 
change: at each step, instead of deleting row i and column j, move them to the far bottom 
and right. The final matrix has determinant Ci ■ ■ ■ e„ det{A); on the other hand, it is diagonal 
with entries congruent to ai, . . . , a„ modulo tt*". 

To obtain the desired estimate, expand det{A) as a sum of signed products ±y4io-(i) ■ ■ ■ Ana{n) 
with a a permutation. When measuring the valuation of the product, each term Aij with 

1 = j contributes f (cj) and each term Aij with i j contributes m; hence for a different 
from the identity, this product has valuation at least the right side of fll.6.3.ip . The product 
All ■ ■ ■ Ann — di^i ■ ■ ■ o-n^n IS the sum of the expressions 

aiCi ■ ■ ■ aj^iCj^iAjj ■ ■ ■ Ann ~ ctiCi ■ ■ ■ ajejA^^jj^i-^i^jj^i-^ ■ ■ ■ Ann 

for j = 1, . . . , n, which have respective valuations m — v{aj) + XliLi '^i^i)- This yields the 
claim. □ 

Remark 1.6.4. Note that Proposition 11.6.31 may be used to obtain approximate character- 
istic polynomials, by applying it with the field K replaced by the completion of K{t) for 
the Gauss valuation (i.e., the valuation which on a polynomial returns the minimum valu- 
ation of any coefficient) and approximating det{tl — A). If one knows that the matrix A 
is "nearly" divisible by some vr*, then one may obtain better information by approximating 
det(7r*t/ — A); this often happens in the setting of p-adic cohomology. For example, suppose 
X is a smooth proper variety over Z^, where g is a power of a prime p < dim(X]F^). Then 
if one writes down the Frobenius matrix on the j'-th rigid cohomology with respect to an 
appropriate basis (namely, a basis of crystalline cohomology modulo torsion), the Hodge 
numbers 

h'^^-' = dimQ, H^-\Xq^, (^ = 0, . . . , j) 

give the multiplicities of as elementary divisors of the matrix. See [231 Theorem 1.3.9] for 
a more general statement. 

Remark 1.6.5. It may be possible to obtain even better bounds on characteristic polyno- 
mials which are more adaptive (i.e., give individual bounds for each coefficient) by using 
more careful linear algebra plus p-adic floating point arithmetic. (In fact, the bound in 
Proposition 11.6.31 is most naturally phrased in terms of floating point arithmetic: the first 
term in (11. 6. 3. II) is the minimum precision of the mantissa of an entry in the final matrix.) 
In particular, it would be interesting to give such bounds for the more general setting where 
the accuracy may vary from entry to entry; in our application to bounding Picard num- 
bers, being able to work in this generality might lend some fiexibility to the cohomological 
calculation. We will not consider the more general setting here. 

2 A little p-adic cohomology 

In this part, we set up a bit of the theory of p-adic cohomology for use later on; this involves 
some consideration of algebraic de Rham cohomology. Some of the calculations, particularly 
Theorem 12.2.51 may be of independent interest. 
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2.1 A homological calculation 

We start with a brief excursion into homological algebra, following [121 Chapter 17]. 

Definition 2.1.1. Let i? be a ring and choose x = [xi, . . . , x„) G -R". The Koszul complex 
K{x) is the exterior algebra A|j(-R") with differentials 

dx{z) = X A z. 

Let K'{x) denote the dual complex, whose underlying i?-module we also identify with 
A^(i?"); let denote the differentials in K'{x). 

Lemma 2.1.2. Let R be a ring. For any x,y & and any z G A|j(-R"), 

{d^dy + dydx){z) = {xiyi H h 

Proof. An easy calculation: see [121 Lemma 17.13]. □ 

Proposition 2.1.3. Let R be a ring, and let C be a complex in the category of R-modules. 
Then for (xi, . . . , x„) G -R", the homology of the product complex K{x) ®C is annihilated by 
the ideal (xi, . . . , x„). 

Proof. Given r in the ideal (xi, . . . , x„), choose ...,?/„ G i? such that xiyi + - ■ ■+Xnyn = f- 
Then Lemma [2.1.21 shows that multiplication by r is homotopic to zero on -ft'(x), with the 
homotopy given by dy. Tensoring that homotopy with the identity map on C yields the same 
assertion on K{x) ® C, proving the claim. □ 

2.2 Algebraic de Rham cohomology 

Algebraic de Rham cohomology is usually considered over a field of characteristic zero, 
but for p-adic cohomological calculations, we also need to work with it over arithmetically 
interesting base schemes. 

Definition 2.2.1. By a smooth pair (resp. smooth proper pair) of relative dimension n over 
a scheme S", we mean a pair (X, Z) in which X is a smooth (resp. smooth proper) S'-scheme 
of relative dimension n, and Z is a relative reduced normal crossings divisor on X. That 
is (in the smooth-only case), etale locally on X, we can find an ^-isomorphism of X with 
a Zariski open subset of under which Z is carried to an open subset of a union of some 
(or all, or none) of the coordinate hyperplanes. We think of Z as defining a logarithmic 
structure on X in the sense of Kato [25], and write (X, Z) also for the resulting log scheme. 
If Z = 0, we abbreviate (X, Z) to simply X. 

Definition 2.2.2. Let (X, Z) be a smooth pair of relative dimension n over a scheme S. 
Let ^{x,z)/s be the sheaf of differentials on X with logarithmic singularities along Z (i.e., 
one allows dt/t for t a local parameter of a component of Z), relative to S; then 
is a locally free coherent Ox-module of rank n. Put ^\xz)is ~ ^hx^(^^^)/^'^• ^^en exterior 
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differentiation induces maps di : Q''^^ — > ^^(^^z)/s ^^^^ that dj+i o rfj = 0. The resulting 
complex is called the de Rham complex of (X, Z) over S", and its j-th hypercohomology 

,Z)IS) 

is called the j-th algebraic de Rham cohomology group of (X, Z) over 5*. 

Remark 2.2.3. Since the f2* are quasi- coherent, we may calculate algebraic de Rham co- 
homology on the etale site instead of the Zariski site [371 Proposition 3.7]. This permits the 
use of etale localization arguments. 

Let (X, Z) be a smooth pair over C, and put U = X \ Z; one then has an isomorphism 

Hi^{{X,Z)/C)^Hi^{U). (2.2.3.1) 

Namely, by Serre's GAGA theorem [12] on the left side and Grothendieck's comparison the- 
orem [20] on the right side (which also uses GAGA, together with resolution of singularities), 
this may be checked for analytic de Rham cohomology, where it amounts to the Poincare 
lemma (see [10] for variations on this theme). 

However, it was pointed out to us by Johan de Jong that one can also establish (12.2.3.11) 
algebraically. In so doing, one can also prove an integral variant where one compares coho- 
mology of the complex of differentials with logarithmic poles with the complex of differentials 
where the poles are made somewhat worse. Here is the result; in its relevance to computing 
p-adic cohomology, it should be viewed as a generalization of [26t Lemma 2]. 

Definition 2.2.4. Let (X, Z) be a smooth pair over a scheme S, and put U = X \ Z. For 

m a nonnegative integer, write (as usual) fi-^^ ^y<^{mZ) for the twist fi-^^ ®Ox ^x{ttT'Z)] 

note that ^{j/g is the direct limit of the Vt'^^ ^y ^{mZ) as m increases. 

Theorem 2.2.5. Let (X, Z) he a smooth pair over a scheme S. For each nonnegative integer 
m, the cokernels of the maps on homology sheaves induced by the natural map of complexes 
of sheaves 



^'{X,Z)/S ^'{X,Z)/si^^) 



are killed by lcm(l, . . . , m) 



Proof. The claim may be verified stalkwise on S, so we may assume S = Spec A is afiine 
and local. We may also localize on X; starting with a point on X, we can shrink to ensure 
that X = Spec R is afiine, ^x/s generated freely by dxi, . . . , dxn for some Xi, . . . ,Xn & R, 
Z = V{xi ■ ■ ■ Xh), and our chosen point lies on V{xi, . . . , Xh)- For each T C {1, . . . , h}, let 
It be the ideal generated by Xi for each i E T, and put Rt = R/It- By a further Zariski 
localization, we may also assume that R contains a copy of each Rx. 
Define 

l<^</i fa;. 9 l<t<h 

dxi = < ^f' . , di= < , 

dXi I > h \ I > h. 
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For each subset U = {ii, ... ,ir} of {1, ... ,n} with ii < ■ • ■ < v, put 

dxu = dxi^ A ■ ■ ■ A dxi^. 
Let M be the set of monomials ^-L tic ^ • • • 

with ji, . . . ,jh e {0, . . . , m}, viewed as 
a partially ordered set under divisibility. For D a nonempty down-closed subset of M (i.e., 
one in which inclusion of fi implies inclusion of any divisor of /i, so that in particular 1 G -D), 
define 

Qd=\J f^-'R; 

note that each di sends Qd into itself, so that Qd ®r z)/s again a complex. 

Let Z) be a down-closed subset of M strictly bigger than {1}. Choose /i = ■ ■ ■ x-'^ 
maximal in D; then D' = D \ {/i} is also down-closed. Let T be the set of i G {1, . . . , /i} for 
which ji 7^ 0, which is necessarily nonempty since we cannot have /i = 1; then the chosen 
inclusion Rx ^ R induces an isomorphism 

/i ^Rt — Qd/ Qd'- 

This isomorphism is equivariant for the action of the di for i ^ T; for i G T, it converts mul- 
tiplication by —ji on the left side into the induced action of di on the right. It consequently 
induces an isomorphism of complexes 

Z{ji ■.ieT)0z ^R^/A - {Qd/Qd') ®r ^{x,z)/s 

= {Qd ®ij ^lx,z)/s)/ (Qd' ®/? ^lx,z)/s)- 

By Proposition 12.1.31 the homology of these complexes is killed by gcd(ji : i 

Now suppose u G Qd ®r ^\xz)/s sa^tisfies duj = 0. Write u = J2u dudxu, where U 
runs over r-element subsets of {l,...,n}. Let hjj G fi^^Rx be the image of gu under 
the map Qd/Qd' 1^~^Rt] then di^jjhudxu) = as well. (It may clarify to think 
of hu as a "leading term" in an expansion of gu as a formal Laurent series in the Xi for 
i G T.) By the previous paragraph, '^^hudxu times gcd(jj : i G T) is a coboundary. Since 
1 < gcd(ji : i G T) < m, C(j is thus equal to a differential whose image in the cokernel of the 
map on homology is killed by lcm(l, . . . ,m) plus a cocycle in Qd' ®r ^\x z)/s- "^^^^ yields 
the claim by induction on the cardinality of -D. □ 

Corollary 2.2.6. Let (X, Z) he a smooth pair over a field K of characteristic 0, and put 
U = X \ Z . Then for each i, the map i7^j^((X, Z)/ K) H\j^{U / K) is an isomorphism of 
K-vector spaces. 

Proof. Apply Theorem 12.2.51 and note that formation of cohomology commutes with direct 
limits. □ 

Remark 2.2.7. Note that if Z is smooth, then the proof of Theorem 12.2.51 actually gives 
something slightly stronger: any closed r-form with poles of order m + 1 can locally be 
written as a closed logarithmic r-form plus a form which when multiplied by lcm(l, . . . , m) 
becomes the differential of an (r — l)-form with poles of order < m. This is crucial for our 
application: see Proposition 13.4.61 
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Proposition 2.2.8. Let {X,Z) be a smooth pair over a scheme S, with Z also smooth; let 
j denote the inclusion Z •—^ X . Then there is an exact sequence of complexes of coherent 
sheaves on X : 

yielding an exact sequence in cohomology 

■ ■ ■ H^^iX/S) -> //^j,((X, Z)/S) ^ H^,^\Z/S) ^ Hli\X/S) ^ ■ ■ ■ . 

Proof. Locally on X we can choose a coordinate x„ cutting out Z; then each section of the 
quotient ^\^x z)/s/^x/s admits a representative of the form dxn/xn A uj for some uo, and 
the residue map carries this section to the reduction of uo modulo Xn- To see that this is 
well-defined globally, we must simply observe that the map is not changed by changing the 
choice of the parameter a;„: if m is a unit on X, then d{xnu) / {xnu) = dxn/xn + du/u, so 
dxn/xn A UJ and d{xnu) / (xnu) A UJ represent the same element of the quotient. □ 



2.3 j9-adic cohomology in general 

We suggest [24j as a useful survey for the material underlying this subsection. 

Notation 2.3.1. Throughout this subsection, let A; be a perfect field of characteristic p, and 
write K for FYa.cW{k) and Ok for W{k). Let a : Ok ^ Ok denote the Witt vector Frobenius, 
which is the unique lift to o^: of the absolute Frobenius on k. 

Definition 2.3.2. Let {X,Z) be a smooth proper pair over k. Let Hl^^^{X,Z) denote 
the i-th (log-) crystalline cohomology of {X,Z); this is an Oj^-module whose construction is 
contravariantly functorial in the pair {X,Z). Moreover, the absolute Frobenius on k acts 
(T-semilinearly on H^^^^^X, Z). For the general construction, see for instance jlSj Chapter 2]. 

To make the Frobenius action on de Rham cohomology explicit, we need to pass to rigid 
cohomology, so we can use the Monsky-Washnitzer interpretation of p-adic cohomology. 

Definition 2.3.3. Let X be a smooth fc-scheme. Let Hl-^J^X) denote the i-th rigid coho- 
mology group of X in the sense of Berthelot [3]; it is a finite-dimensional i^- vector space 
whose construction is contravariantly functorial in X. Moreover, the absolute Frobenius on 
k acts cr-semihnearly on if*ig(X). For X proper, Berthelot [U Proposition 1.9] constructs a 
functorial, Frobenius-equivariant isomorphism 

For (X, Z) a smooth pair over fc, Shiho ^Hl §2.4] (using crucially a result of Baldassarri and 
Chiarellotto [2]) constructs a functorial, Frobenius-equivariant isomorphism 

Hl(X,Z)®,^K^HUX\Z). 
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2.4 j9-adic cohomology in the liftable case 

When things can be hfted nicely to characteristic zero, the construction of ]9-adic cohomology 
becomes much simpler. 

Proposition 2.4.1. Let (X, 3) be a smooth proper pair over Ok- Put X = X^, Z = 3fc; O'lT'd 

u = x\z. 

(a) There are isomorphisms 

Hl,^,iX,Z) = W^^iiX,3)/0K) 

which are functorial in smooth proper pairs over Ok- In particular, the right side 
inherits an action of Frohenius. 

(h) Suppose that 3 is also smooth, and use Proposition \2. 2. 8\. Definition \2.3.3\. and (a) to 
obtain an exact sequence 

> H^ig{X) H^ig{U) Hlr^{Z) if*+^(X) ^ • • • . 

These maps are then Frobenius-equivariant for the following twists: 

Hl,^{X) ^ Hl,^{U) 
Hl,^{U) ^ Hlr^\Z){-l) 
Hl^\Z)^Hlt,\X){l), 

where M{n) denotes M with its absolute Frohenius action multiplied hy p~^ . 

Proof. For (a), see [25l Theorem 6.4]. For (b), we may invoke rigid analytic GAGA to argue 
that algebraic de Rham cohomology of a smooth proper pair over K coincides with rigid 
analytic de Rham cohomology of the analytification. The Frohenius equivariance may now 
be checked at the level of complexes by following the construction of the Gysin isomorphism 
in rigid cohomology P, §5], [16]. □ 

One also has a nice description of the Frohenius action on the rigid cohomology of a 
smooth affine scheme in terms of a lifting. 

Definition 2.4.2. Let X be a smooth affine A;-scheme, suppose X is a smooth affine Ok- 
scheme lifting X, write X = Spec A, and choose a presentation A = Ok[xi, . . . Let 
Ok{xi, . . . , XnY be the ring of power series in Oii:|xi, . . . , which converge on some polydisc 
of radius greater than 1, and put 

A"^ = Ok{Xi, . . . ,Xny /IOk{Xi, . . . ,X„)^ 

Then Berthelot [3^ Proposition 1.10] constructs an isomorphism between Hl-^J^X) and the 
cohomology of the complex 

the latter is the complex computing Monsky-Washnitzer's "formal cohomology" [17]. More- 
over, a F : A'^ ^ A^ is any ring homomorphism extending a on Ok and lifting the absolute 
Frohenius, then pullback by F induces the Frohenius action on rigid cohomology. 
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Remark 2.4.3. For example, one can compute from the above description that the g-power 
Frobenius action on H'^^^{F^J consists of multiphcation by q\ This also follows from the fact 
that the cohomology of projective space is generated by cycle clasess. 



3 The case of smooth hypersurfaces 

From the first part of this paper, we obtain a procedure for bounding from above the Picard 
number and the geometric Picard number of a smooth proper variety X over a finite field 
Fg-. compute an approximation modulo p"^, for some m, to the matrix via which Frobenius 
acts on the rigid cohomology space H'^{X) over Qg (i.e., a "higher Cartier matrix"), then use 
Algorithm 11.6.21 to bound the right-hand side of fll.5.2.ip or f ll.5.2.2p , respectively. What 
remains to be done, for any particular class of varieties, is to describe how to compute 
the approximate Frobenius matrix for varieties in that class, by realizing the constructions 
described in the second part of the paper. Here, we describe one such procedure for smooth 
hypersurfaces in a projective space, based on work of Griffiths [TT], and give a few details 
of an implementation of this procedure which we have constructed. We also mention some 
alternate approaches that we have not (yet) experimented with. 

It is worth noting that much of what is described below generalizes relatively easily to 
smooth hypersurfaces, or even smooth complete intersections, in toric varieties; for example, 
Johan de Jong is currently developing an implementation for hypersurfaces in weighted 
projective spaces [8]. We have restricted to hypersurfaces in projective space, and ultimately 
to surfaces in P^, both for simphcity of exposition and because that is all that we have 
attempted to implement ourselves. 

3.1 A calculation on projective space 

We pause for a brief excursion into the cohomology of sheaves of differentials on projective 
space. 

Lemma 3.1.1. Let S be a scheme. For any positive integer n, there is an exact sequence of 
sheaves on X = Pg.- 

Proof. A standard calculation: see [2TI Theorem II. 8. 13]. □ 

The following is due to Bott |6] over C (and can be deduced over any field of charac- 
teristic zero); we were unable to find a reference for the general version, so we include the 
easy calculation. One might like to think of it as a special case of the Kodaira-Nakano 
vanishing theorem, but the latter is not valid over a general base; Raynaud [ID] exhibited a 
counterexample over a field of positive characteristic. 

Proposition 3.1.2. Let S = Spec A be an affine scheme. For n a positive integer, put 
X = Pg. Let k,p,q be integers with < p,q < n. Then H'^{X,Q^^^{k)) = unless one of 
the following conditions holds: 
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(a) p = q and k = 0, in which case H^{X, f2^y^(A;)) is free of rank 1 over A; 
(h) q = and k > p; 
(c) q = n and k < p — n. 

Proof. We will use without comment Serre's calculation of the cohomology of Ox{k), in the 
form of the statement of the lemma in all cases where p = 0; for this, see Theorem III. 5.1] 
(and suppress the superfluous noetherian hypothesis) or Proposition 2.1.12]. 

Take the p-th exterior power of the exact sequence in Lemma I3.1.H then twist by k, to 
obtain 

- n^^^,{k) {A^{Ox{-ir^')) (k) - n^^-^lik) 0, 

which we may rewrite as 

^ ^x/sik) Ox{-p + k)^"p''> ^ ^""x/lW 0- (3.1.2.1) 
We now proceed by inspecting part of the long exact sequence in cohomology of fl3.1.2.ip : 

H^-\X,Ox{-p+k)ft') ^ /J^-i(X,(]^-^(fc)) ^ iJ''(X,fi^/^(A;)) ^ H^iX,Oxi-p+k)ft') 

(3.1.2.2) 

In case k < p—q and < q < n, then the outside terms of fl3.1.2.2p vanish, so H'^~'^{X, f2^^(fc)) 
and H'^{X, Q^^^^^k)) are isomorphic. We obtain the vanishing of H''{X, fl^^^^^k)) in this case 
by successively decreasing both p and q until the step when at least one of them goes negative, 
at which moment the correct formal interpretation of f l3.1.2.2p yields the desired vanishing. 
Similarly, in case k > p — q and < g < n, we obtain vanishing of H'^{X,Q^j^^g{k)) by 
successively increasing both p and q until the step when at least one of them exceeds n. 

In case k = p — q, the outside terms of (13.1.2.20 still vanish as long as 1 < g < n. If 
k 0, then successively decreasing p and q still eventually manages to send one of them 
below zero, so we get the desired vanishing. If = 0, we instead hit H^{X,Ox) which is 
free of rank 1. This yields all of the desired results. □ 

Corollary 3.1.3. Let S = Spec A be an affine scheme, let n be a positive integer, and put 
X = Pg. Let k be a nonnegative integer, let Z be a smooth hypersurface in X , and define the 
complex = ^^^^((fc + p)Z) with the evident differentials d. Then the hypercohomology of 
C' coincides with the homology of the complex H^{X,C'). In particular, 

H"(X,C-) = if°(X,C")/c/if°(X,C""i). 

Proof. We compute EI"(X, C ) using a spectral sequence with Ef^ = if^(X, C^). By Propo- 
sition EU and the fact that ^ ^]^/^((A; + |j) deg(Z)), we have = for g > 0. Hence 
the sequence degenerates at E2 and yields the claim. □ 



Corollary 3.1.4. For any affine scheme S = Spec A, 



A i = 0,2,...,2n 
otherwise. 



16 



Proof. By Proposition 13.1 .21 H^{F^, f2p„y^) is free of rank 1 over Aii p = q & {0, . . . , n} and 
is otherwise zero. Hence the spectral sequence computing hypercohomology degenerates at 



3.2 Cohomology of smooth hypersurfaces (after Griffiths) 

The middle cohomology of a smooth hypersurface in a projective space was described by 
Griffiths [T7] using mostly analytic arguments (i.e., working over C and invoking GAGA). 
One can reconstruct these results algebraically; we will not do so explicitly, but we will use 
algebraic techniques later to extract arithmetic information. 

Notation 3.2.1. Throughout this section, let be a field of characteristic zero, put X = 
P^, let Z he a smooth hypersurface of degree din X, defined by the homogeneous polynomial 
P G K[xo, . . . ,Xn], and put U = X \ Z . By the Lefschetz hyperplane theorem, the map 
H^j^{X/K) — » H\^{Z/K) induced by the inclusion Z X is bijective for i < n — 2 and 
injective for i = n — 1. By Corollary 13.1.41 and Poincare duality, we thus conclude that for 
< i <2n — 2 with i ^ n — 1, we have 



In particular, the only cohomology group of Z requiring further inspection is the middle 
cohomology H2^\Z/K). 

The following result lets us shift attention from Z to the ambient projective space X, 
where it is easier to make calculations. Remember that ifjf^((X, Z)//^) = H^Yi{U/K) by 
Corollary [2Xni 

Proposition 3.2.2. The exact sequence of Proposition \2.2.S\ induces exact sequences as 
follows: if n is even, then 



El, yielding the desired result. 



□ 




^ H^^{U/K) H2^\Z/K) 



is exact; if n is odd, then 



HUU/K) ^ H2^\Z/K) H2l\X/K) 



is exact. 



Proof. See [HI (10.16)]. 



□ 



Definition 3.2.3. Put 



n 




{—lyxi dxo A ■ ■ ■ A dxi A ■ ■ ■ A (ix 
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where the hat denotes omission. One then calculates as in \T7, §4] that H2ji{{X , Z) / K) is 
isomorphic to the quotient of the group of n-forms AQ/P"^, where m is a positive integer 
and A E K[xo, . . . , Xn] is homogeneous of degree md — n — 1, by the subgroup generated by 

(diA)n A(diP)n 
— 

p-m pm+1 

for all nonnegative integers m, all i G {0, . . . , + 1}, and all homogeneous polynomials 
A G K[xq, . . . , Xri\ of degree md — n. (Here di is shorthand for ^.) 

Remark 3.2.4. Note that H^^iU / K) admits a natural filtration whose i-th step consists of 
those classes represented by forms Afi/P™ for some integer m <i + l. In fact, this filtration 
is induced by the Hodge filtration on H2^\Z/K) [TT, §10]. 

Remark 3.2.5. The description in Definition 13.2.31 gives rise to a natural "reduction of 
poles" procedure for computing in H^^{U), sometimes referred to as the Grifiiths-Dwork 
method. First, one writes down a basis: for h = 1, . . . ,n, one finds monomials of degree 
hd — n — 1 which generate the quotient of the space of all such monomials by the multiples 
of doP, . . . , dnP- Then, to write the class of a given form AQ/P"^ in terms of these, one uses 
a Grobner basis division procedure to write A as a linear combination of OqP, . . . , dnP (plus 
basis elements if m < n), then reduces the pole order. The fact that it is always possible to 
perform this reduction follows from a theorem of Macaulay pT| §4] or from a sheaf-theoretic 
reinterpretation [17, §10]. 

3.3 The p-adic cohomology interpretation 

We now use the previous subsection to describe the p-adic cohomology of a smooth hypersur- 
face in F^^. Note that while the p-power Frobenius action on rigid cohomology of a variety 
over Fq will only be semilinear, the g-power Frobenius action will be linear over Qq. 
We start by setting notation for the rest of the section. 

Notation 3.3.1. Let Z he a smooth hypersurface of degree d in X = F^^, for Fg a fi- 
nite field of characteristic p > 0, defined by the homogeneous polynomial P{xq, . . . ,Xn) G 
Fq[xo, . . . , Xn]- Choose a lift ^(xq, . . . , x„) G ^^[xo, . . . , Xn] of P to a homogeneous polyno- 
mial of the same degree d, let 3 be the zero locus of ^ in X = P^^, and put ii = X \ 3- Put 
X = Xq^, Z = 3q,, and U = iiq^; also write P for ^ as an element of Qq[xo, . . . , x„]. 

By Definition [2X3l H^igiZ) = Hl^^{Z) K- adding Proposition EXIKa), we get 
Hi.^{Z) ^ i/*j^(3/Zg) Qg = H\^{Z/qq). Adding the Lefschetz hyperplane theorem, 
we obtain Hl^^{Z) = if*ig(X) for i < n — 2. Adding Poincare duality, we see that for 
< i < 2n — 2 with i n — 1, if i is odd, then Hl^^{Z) = 0, while if i is even, then Hl-^J^Z) 
is one-dimensional and the g-power Frobenius acts by multiplication by g*/^. Moreover, by 
Proposition 13.2.21 and Proposition 12.4. 1( b) . we have Frobenius-equivariant exact sequences 

Q ^ H-^iU) ^ H:r\Z){-l) ^ Q 
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if n is even and 

^ H^,^iU) ^ H:^-\Z){-1) H-^\X) ^ 

if n is odd. That is, i7^g(f/)(l) coincides with except that if n is odd, its generahzed 

eigenspace for Frobenius of eigenvalue g^""^)/^ has dimension one less. (That is, if^"g(f/) is the 
primitive middle cohomology of Z, i.e., the part not explained by the Lefschetz hyperplane 
theorem.) 

Thanks to Berthelot's comparison theorems (see Definition I2.4.2p . we can describe the 
Frobenius action on H'^^JJJ) as follows. 

Definition 3.3.2. Let v denote the Gauss valuation on polynomials over Q^; that is, 
vi^cix^) = min/{fp(c/)}, where Vp denotes the p-adic valuation on normalized by 
Vp{p) = 1. 

Definition 3.3.3. Let R denote the ring of formal sums Xli^o where Ai G Qq[xo, . . . , x„] 

is homogeneous of degree di, and 

hminf^^^>0. 

That is, the valuations of the grow at least linearly in i. Define the ring map F : R ^ R 
by formally setting F[ tic ^ ^ oc ^ for i = 0, . . . , n, and 

expanding the parenthesized expression as a binomial series. Formally extend F to n-forms 
by setting 

F{An) = F{Axo ■ ■ ■ x„)F(xo 1 ■ ■ ■ x^^fi) , 

where F{dxo/xo) = q dxo/xo and so forth. As noted in Definition l2.4.2[ this ring map induces 
the g-power Frobenius in rigid cohomology on H^,g{U) = H'^^{{X, Z)/%) = H^^{U/Qq). 

3.4 Precision estimates 

The plan now is to compute an approximation to the Frobenius action on H^^^iU) by ap- 
plying a truncation of F to a basis of H^^^iU) and using the "reduction of poles" process 
(Remark I3.2.5P to put the results back in terms of the basis. To do this, we must produce 
effective bounds on the amount of p-adic precision needed to keep the error introduced by 
the truncation to a particular amount. One can derive general bounds easily from the theory 
of p-adic cohomology, but for effective bounds we must work a bit harder. This is analogous 
to the precision analysis in but the higher-dimensional situation we are in makes things 
a bit more technical. 

Our first order of business is to relate a basis, which is natural from the point of view of 
reduction of poles, to the integral structure on de Rham cohomology. 
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Definition 3.4.1. Let H be the image of H^^{{X,3)/^q) in H2J^{{X , Z)/Qgy, we refer to 
elements of H as integral elements of Z)/Qg), or of H'2^{U /Qg), or of H^.giU). 

Definition 3.4.2. Let B denote a basis of H2j^{U /Qg) obtained as follows. For h = 1, . . . ,n, 
form the quotient of the space of homogeneous polynomials in ¥q[xo, . . . , x„] of degree hd — 
— 1 by the multiples of d^P, . . . , dnP- Find monomials m¥g[xo, . . . , x„] which project onto 
a basis of this quotient, then lift these monomials to monomials in Zq[a;o, • • • , Xn]- For each 
such lift /i, include /xfi/P'^ in B. Let V C H^^iU /Qg) denote the Z g-span of B. (This is not 
the only logical choice; see Remark I3.4.12I ) 

Lemma 3.4.3. Let W be the Zg-span in H2^{U /Qg) of the elements hVl/P^', for h E 
{1, . . . ,n} and /i G Zq[xo, • • • , Xn] homogeneous of degree hd — n — 1. Then 

H CW. 



Proof. Map the complex ^'(^x3)/z ^^^'^ complex C' with 

then invoke Corollary 13. 1.3[ □ 

Remark 3.4.4. In general, we have H ^ W hy Lemma [3.4.3^ V evidently, and [n — 
1)\W C V by inspection of the reduction process. In the special case p > n, though, it will 
follow from Corollary 13 . 4. 71 that H = V = W. In this case, we also know that H^^^i^/Zg) is 
torsion-free because of the degeneration of the Hodge-de Rham spectral sequence. See p4] 
for further discussion. 

We next consider the loss of p-adic precision incurred when one reduces a given differential 
into standard form. 



Definition 3.4.5. For m a positive integer, let f{m) be the smallest integer t with the 
following property: for each form u = AQ/P"^ with A G Zg[xo, . . . , x„] homogeneous of 
degree md — n — 1, p^io represents an element of B.. The following relations are evident 
(using the relations in cohomology): 

/(i) = o 

/(m) </(m+l) 

/(m) = f{j>\vn/p)\) [m > n). 

Proposition 3.4.6. For each m > 0, 

n 

f{m) < ^[logpmax{l,m - 

i=l 



20 



Proof. Apply Theorem I2.2.5t we get the n terms of the sum by feeding the cohomol- 
ogy sheaves into the spectral sequence computing hypercohomology, keeping in mind Re- 
mark [2221 (We don't get an (n + l)-st term because the map of Theorem 12.2.51 on zero-th 
cohomology sheaves is an isomorphism, so those do not contribute.) □ 

Corollary 3.4.7. We have f{m) < Vp{{m — 1)\) for allm > 0. In particular, in the notation 
of Lemma\34^ we have (n - 1)\W C H CW. 



This bound is asymptotically nlogp(m), which for our application to surfaces will be a 
bit too large to be practical. Fortunately we can shave it down a bit. 

Definition 3.4.8. For m,i integers with i > 0, let g{m,i) be the p-adic valuation of {~^)- 
By a standard argument attributed to Kummer, g{m,i) equals the number of carries in the 
base p addition of i and — m — i. 

Proposition 3.4.9. Let m be a positive integer. Put 

N = max{/((m + i)p) - i- g{m, £)}. 

Then 

f{mp) < max{A^, n — 1 + f{m)}. 

Proof. Given a form AQ/P"^^ with A G Zq[xo, . . . , a;„], separate the monomials of A depend- 
ing on the reductions modulo p of their exponents. Let B be the sum of the monomials 
whose exponents are all congruent to p — 1 modulo p, and put C = A — B. Then CQ/P"^^ 
is cohomologous to a form mpDQ / P"^^"^^ with D G Zg[xo, . . . , because 

xfn ^ mp xi'^\d,P)n 

pmp J _|_ ]^ pmp+1 

in cohomology. Such a form reduces to an element of p'^H for t = f{mp+ 1) — 1 — g{m, 1) < 
A^. On the other hand, we can write B = p~"F{D) + E with D having pole order m 
and E = ^^>o/(-7)E^fi/p(™+^)p for some E^ G Zjxo,...,xJ. Since p-^F acts on H 
(by comparison with if^j^^(3/Zq) via Proposition 12.4.1]) . D reduces to an element of p~^H 
for t = n — 1 + fim), while the i-th summand of E reduces to an element of p~^H for 
t = f{{m + €)p) — I — g{m, €) < N . This yields the claim. □ 

Proposition 13.4.91 plus any sublinear bound on / suffices to give an upper bound of the 
form (n — 1) logp(m) plus a constant. However, for implementation purposes, it is important 
to control that additive constant as much as possible; this can be done using an iterative 
algorithm as follows. 

Algorithm 3.4.10. Define the function 

n 

Mm) = J^[logpmax{l,m -i}J. 

i=l 
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Given an input positive integer m and an input array Aq such that f{j) < Aq^J) for each j 
for which Ao{j) exists, the following algorithm, if it terminates, returns an array A of length 
at least m, such that f{j) < A{j) for each j for which A{j) exists, and A{j) < Ao(j) for 
each j for which A{j) and Ao{j) both exist. 

1. Create an array 

A( ■\ J^o(j) ^o(j) exists 

Mj) = . (j = l,...,min{n,m}). 

I JoUj otherwise 

2. Make a copy A' of A. Put j = n + 1. 

3. If A{j) is not defined, check whether A and A' are identical (of the same length). If 
so, return A and STOP. Otherwise, go to step 2. (If A{j) is defined, continue to step 
I) 

4. Put ji = p\j/p], N = n — 1 + A{ji/p), and i = 1. 

5. If np < ji + ip and n logp(ji + Ip) — i < N, then set 

Aij) = A{j + l) = --- = A{j,) = mm{N, /o(ji)}, 

replace j by ji + 1, and go to step 3. 

6- If /o(ji + ^p) — (■ — qUiA) ^ N , 'then replace i by i W and go to step 5. 

7. Extend A using the formula A{i) = fo{i) if needed to ensure that A{ji + ip) exists. 
Replace N by max{A{ji + ip) — i — g{ji,i), N}, replace i by i + 1, and go to step 5. 

Proof. What we show is that at every stage, whenever some A{j) is defined, we have f{j) < 
A{j). This holds whenever an A{j) is instantiated by Proposition I3.4.6I In step 5, the 
quantity n logp(ji + ip) — £ is a decreasing function of i for np < ji + ip; it is also an upper 
bound for /(ji + ip) — i — g{ji, i) by Proposition 13.4.61 again. The property that /(j) < A{j) 
is preserved in step 5 thanks to Proposition 13.4.91 □ 

Remark 3.4.11. One can prove termination of Algorithm 13.4.101 and control its runtime 
with a bit of effort, but in practice it suffices to simply let it run until either the process 
terminates, or one goes through a set number of iterations of step 3 (we used 20 iterations 
in our examples). 

Remark 3.4.12. If in Definition 13.4.21 we had taken B to consist of the elements [h — 
ly.fin/P^ instead oifin/P^, we would get a Z^-span V satisfying {n-l)\W CV' C H CW. 
However, in practice this choice appears to give slightly less control of the precision loss. 
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3.5 Summary of the algorithm 

To summarize, we describe how to assemble an algorithm for computing an approximate 
Frobenius matrix on if^"g(f/)(— 1). 

To start with, write down a basis B of rigid cohomology as in Definition 13. 4.21 Say we 
want to compute the Frobenius matrix on this basis modulo p"^. Choose the integer s using 
the following algorithm: 

1. Put s = r. Create an empty array A. 

2. Put j = s-n + 1. 

3. If j > and nlogjj{p{n + — l)<n — — r, then return s and STOP. 

4. Replace A by the result of Algorithm 13.4.101 applied with inputs p{n + j) and A. If 
A{p{n + j)) >n — — r, then replace s by s + 1 and go to step 2. Otherwise, 
replace j by j + 1 and go to step 3. 

For each basis element, compute the image of Frobenius truncating all terms which 
vanish modulo p*, then use reduction of poles (Remark I3.2.5j see also Remark 13.5.21 be- 
low) to express the result in terms of the basis elements. The verification that this gives 
enough precision is straightforward: the Frobenius image of a monomial fifl/P^ has the form 
Y^^=qP^^^~^ BjVt/ P^^^'^^\ so we need to ensure that 

n — 1 + j — f{p{h + j)) > r {h < n, j > s — n + 1) , 

but checking for h = n implies the same for h < n, and the upper bound from Proposi- 
tion [3]3]6] allows us to truncate in step 3. 

Remark 3.5.1. This algorithm readily admits some parallelization, as one can compute the 
reductions of the Frobenius images of different basis elements independently. The experi- 
mental results we describe later do not depend on this capability, but it may be useful for 
larger examples. 

Remark 3.5.2. There are at least four distinct ways to carry out the reduction of poles 
implied by Remark 13. 2. 5[ 

• One can simply perform the entire calculation over Q, then interpret the final result 
modulo the appropriate power of p. This was our first choice, implemented in Singu- 
lar, but it leads to undesirable intermediate coefficient blowup. 

• One can perform the calculation over Z/p^Z by using integral analogues of Grobner 
basis arithmetic as implemented in Magma: to do this, one must postpone the division 
by m — 1 implied in reducing the pole order from m to m — 1 until the end of the 
calculation. One must also remember to use the trivial relation AQ/P"^ = APQ / P"^^^ , 
as P may not be generated by its partial derivatives (in case p divides deg(P)). This 
was our second choice, and is used in our current implementation. 
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• One can perform the calculation over for any prime ^ of good reduction of the lifted 
hypersurface. This would allow the use of the more efficient polynomial arithmetic of 
Singular over Magma; however, in order to recover the desired answer (by working 
modulo many small i and using the Chinese remainder theorem) one would need a 
bound on the heights of the entries of the resulting matrix. We have not attempted 
this method. 

• One can perform the calculation over C by going through the Betti-de Rham compar- 
ison as in p/fj: numerically integrate each truncated Frobenius image against a basis 
of homology, then perform a lattice reduction to express these in terms of the integrals 
of a basis of cohomology. Again, this requires height bounds, and again we have not 
attempted this method. 

The last two methods share ideas with the method used by Edixhoven et al [TT] to compute 
coefficients of the A modular form (and by extension with the Schoof-Elkies-Atkin method 
for computing zeta functions of elliptic curves). 

3.6 Alternate algorithms 

As the experimental data in the final part of the paper suggests, computing approximate 
Frobenius matrices in p-adic cohomology in the manner we have suggested is rather laborious. 
There are several other ways one might perform this computation; we mention some of these 
in passing, noting that any of them can be used together with the first part of the paper to 
give a test for low Picard number. (The arguments in the second part of the paper, notably 
Theorem 12.2.51 may help in the analysis of precision loss in such algorithms.) 

The shift from a hypersurface to its affine complement amounts to an increase by one in 
the dimensions of the varieties under consideration, and in the number of variables in the 
polynomial rings in which the calculations take place. (In reality there is one more variable 
even than that, but this is merely because we are working with homogeneous polynomials.) 
That shift turns out to be costly, so one would ideally like to avoid it. This appears to be 
possible for so-called nondegenerate hypersurfaces, those which together with the coordinate 
hyperplanes and the hyperplane at infinity form a normal crossings divisor. For curves, 
this has been proposed by Castryck, Denef, and Vercauteren [7], but the method extends 
readily to higher dimensions. We made a cursory attempt to implement this method for 
surfaces in P^, but our results were inconclusive: the additional complexity in the method 
(especially in lifting Frobenius, which would now be done on an affine piece of the original 
hypersurface rather than on an affine complement) seemed to introduce large constant factors 
that interfered with the asymptotic improvements at the scale at which our calculations took 
place. Nonetheless, we think the method deserves further study. 

A better approach may be to use devissage: write the given surface as a fibration of curves, 
compute the higher direct images of the constant sheaf, then compute the cohomology of 
these. This has been suggested by Lauder, who has implemented this in some examples with 
good results [32] . 
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Yet another approach is to avoid directly computing the cohomology of the particular 
hypersurface of interest, by instead putting it into a pencil with one member chosen to be 
smooth with extra automorphisms. One can then compute its Frobenius matrix more easily, 
then use that data as the "initial condition" in the differential equation relating the Gauss- 
Manin connection of the pencil with its global Frobenius action. This is the "deformation" 
method of Lauder [301 El]; it has been tested experimentally for families of elliptic curves 
by Gerkmann [15] and has been theoretically analyzed for hyperelliptic curves by Hubrechts 
[23] (where it already gives some improvement over the direct method), but we are not aware 
of any work in higher dimensions besides Lauder's original papers. 

4 Implementation details 

In this section, we describe an implementation that implements a special case of the algorithm 
we have described, and give some experimental results. One glaring omission is that we do 
not make a complexity analysis; this is only partly out of laziness. The other reason is 
that Grobner basis calculations in general have extremely bad worst-case performance; we 
are not in the worst case here, but we would have to look closely at what we are using to 
obtain complexity estimates. Since the purpose of this paper is instead to demonstrate the 
practicality of these methods, we do not carry out such intricate analysis here. 

4.1 Implementation notes 

Using the Magma algebra system [S], we have developed an implementation of the methods 
of this paper, to obtain an algorithm for computing approximate Frobenius matrices in rigid 
cohomology for smooth surfaces in projective 3-space over a prime finite field. See [1] for the 
source code. 

We have tested this implementation on the computer dwork . mit . edu, a Sun workstation 
with dual Opteron 246 CPUs running at 2 GHz, with access to 2GB of RAM. Although 
these CPUs are 64-bit processors, these experiments were conducted in 32-bit mode under 
Red Hat Enterprise Linux 4. Each individual surface was run on a single CPU with no use 
of parallelism (see Remark 13.5. II) . and timings are reported in CPU seconds; memory usages 
are reported in megabytes. Beware that these should only be taken as order-of-magnitude 
indications: there are slight variations from run to run of a single example, and there are 
much bigger variations within classes of examples, probably arising from the vagaries of 
Grobner basis arithmetic. 

Some initial experiments were also conducted using the Singular algebra system |16j . 
The main downside with using Singular for this calculation is that it only treats poly- 
nomials over a field; while one obtains correct answers by reducing poles over Q and then 
reducing modulo a power of p, the resulting calculations experience unacceptable interme- 
diate coefficient blowup. By contrast, the Magma implementation uses a Grobner basis 
implementation over Z/p^Z, which avoids the coefficient blowup. (Compare Remark 13.5.21 ) 
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In the subsequent sections, we describe some examples computed using this implemen- 
tation. These examples were chosen to be "generic", without special geometric properties: 
their coefficients were chosen at random with a bias towards zero coefficients. The bias 
somewhat simplifies the Grobner basis calculations. 

It is worth noting that one can use a "prescreening" strategy to find such examples: 
deliberately compute approximations with not enough initial precision, then revisit the ex- 
amples that appear to work with a provably sufficient amount of precision. We suspect that 
this works because our precision estimates do not give a complete picture of how quickly the 
p-adic approximations are converging; see Remark 14.2.31 for an instance of this. 

4.2 Example: degree 4 over F3 

We start with a careful analysis of an example in what is possibly the simplest nontrivial 
case. Namely, surfaces of degree 1 and 2 over a finite field Fp are isomorphic to and 
X P^, whose zeta functions are known, while surfaces of degree 3 have all cohomology 
generated by the classes of the 27 lines on the surface over Fp, so the zeta function can be 
computed from the Galois action on these lines. (Note that while a smooth cubic surface 
over Fp has geometric Picard number 7, its arithmetic Picard number can equal 1; see |53j 
for an example.) 

A smooth surface in P^ of degree 4 is a K3 surface, whose middle cohomology has dimen- 
sion 22 and Hodge numbers 1, 20, 1. (Remember that we compute using primitive middle 
cohomology, in which the dimension and the central Hodge number are both decreased by 
1.) We will exhibit an example of provable arithmetic Picard number 1 over F3; the more 
natural first choice F2 actually turns out to be somewhat trickier (see Section |4^ . 

Example 4.2.1. The smooth quartic surface over F3 defined by the polynomial 

— xy^ + xy^w + xyzw + xyw"^ — xzw'^ + y^ + y^w — y'^zw + z^ + w'^ 

was found to have Picard number 1 by computing a Frobenius matrix modulo 3'^. To carry 
out this calculation provably using the optimal bound extracted from Algorithm 13.4. 10[ it 
was necessary to truncate differentials modulo 3''; using only Proposition 13.4.61 would have 
required working modulo 3^^. 

Remark 4.2.2. For comparison. Table [T] gives some timings and memory usages for various 
levels of initial and final precision in Example 14.2. 1[ When a final precision is specified, 
that means the initial precision in that row is the minimum needed to guarantee that final 
precision in the Frobenius matrix under Algorithm 13.4. 10( however, see Remark 14.2.31 

Remark 4.2.3. One reality check on Example 14.2. II is to verify the implied compatibilities 
between the answers computed to various p-adic precisions; that is, the computed Frobenius 
matrices for any two rows in Tabled] should agree modulo the final precision of the earlier row. 
This turns out to be true in a strong fashion: some of the calculations are even more accurate 
than predicted. Namely, we obtain the correct Frobenius matrices modulo 3^, 3^, 3^, 3^ using 
initial precisions 3^, 3^, 3^, 3^°, respectively. 



26 



Table 1: Timings for a smooth quadric over F3. 



Final precision 


Initial precision 


CPU sec 


MB 




36 


227 


37 


33 


37 


731 


53 




38 


907 


64 




39 


4705 


124 


34 


310 


13844 


906 


35 


3" 


15040 


1103 


36 


312 


40144 


1795 



Remark 4.2.4. Another reality check on Example 14.2. H is a comparison of initial coefficients 
of the predicted zeta functions against those coefficients determined exactly by actually 
enumerating rational points. Over W^i for z = 1, 2, 3, 4, 5, the respective numbers of rational 
points on the surface are 

8,80,713,6836,58868; 

this counting took several hours on a laptop computer using a simple-minded MAGMA pro- 
gram. From the formula for the zeta function (Definition 11.4. ip . and the fact that it has 
the form {Q{T){1 - T)(l - 3T)(1 - 9T))-^ where Q(T) has degree 21, we deduce that the 
characteristic polynomial of 3^^F on primitive middle cohomology begins 

-(3r2i + 5r2o + 6Ti9 + 7ri^ + 5T^'' + 4ri^ + ■■■)• 

3 

On the other hand, applying Remark 11.6.41 to the Frobenius matrix computed with final 
precision 3^, we determine that the same characteristic polynomial is congruent modulo 3^ 
to 

-(3T2i + 5T20 + 6T^^ + 7T^^ + 5T^^ + 4T^^ + 2T^^ - T^^ - 3T^^ - 5T^^ - 5T^^ 
3 

-5r^° - 5T^ - 3T^ - + 2r^ + 4T^ + 5T^ + 7T'^ + 6T^ + 5T + 3). 

Not only are the two assertions consistent, but the characteristic polynomial computed from 
p-adic cohomology demonstrates the symmetry forced by the functional equation of the zeta 
function, i.e., by Poincare duality on cohomology. (The geometric Picard number in this 
instance appears to be 4, as the characteristic polynomial appears to be divisible not only 
by T + 1 but also by -|- 1, so Tate's conjecture predicts extra cycle classes over F32 and 
again over F34.) 

Remark 4.2.5. It seems an interesting question to explore to what extent one can recover 
a zeta function from the various pieces of data we have in the above situation: 

• the symmetry and location of roots, from the Weil conjectures; 
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• the initial point counts; 

• divisibilities implied by the relationship between the Newton and Hodge polygons (see 
Remark [LMD; 

• congruences derived from computing p-adic cohomology to low precision. 

In particular, in many cases it may be possible to combine information to recover zeta func- 
tions using p-adic cohomology calculations at much less precision than would be predicted 
by a straightforward application of the Weil conjectures plus taking into account the Hodge 
polygon. In the case of Example 14.2.11 this is discussed in detail in [28] ; the end result is 
that not only is the zeta function equal to the guess of Remark |4.2.4[ but in fact this already 
follows from the computation of the characteristic polynomial of 3~^F modulo 3^ (and so 
from the Frobenius matrix with final precision 3^) without counting any rational points at 
all. 

4.3 Examples: degree 4 over ¥p = 5, . . . , 19) 

We next exhibit some examples where we can compute the geometric Picard number. These 
examples appear in a construction of van Luijk [48, Proposition 5.1]; for context, we first 
state [m Proposition 4.1]. 

Proposition 4.3.1 (van Luijk). Let k be a field, and suppose a,f3 & k satisfy a^(3 ^ a(3^. 
Let f G k[w, X, y, z] be a homogeneous polynomial of degree 3, such that either the coefficients 
ofy^ and in f differ, or the coefficients ofx^y andx^z in f differ. Suppose that the surface 
X in defined by 

wf — {xy + xz + ayz){xy + xz + [3yz) (4.3.1.1) 

is smooth with geometric Picard number 2, and put X = X Xj^ k. Then the group Aut(X) is 
trivial. 

Note that the surface in f l4.3.1.ip . if smooth, has arithmetic Picard number at least 2, as 
the hyperplane section w = splits into two conies. 

Example 4.3.2. Over F5, we verify the instance of Proposition 14.3.1) with a = 1, /? = 3, 
and 

/ = 3x + 3xy — xyw + 3xzw — xw + y — y w + 2z + w ; 

one case of ^E, Proposition 5.1] relies on the fact that this surface has geometric Picard 
number at most 2 (and hence has Picard number and geometric Picard number exactly 2). 
This we verify in turn by computing the matrix modulo 5^, which requires initial precision 
5^, 5482 CPU seconds, and 514 MB of memory; applying Algorithm 11.6.21 to bound the 
contributions to the right side of (11. 5. 2. 21) shows that the geometric Picard number is indeed 
at most 2. 
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For p > 5, there are relatively few values of n that can contribute to fll.5.2.2p for which 
Cn is p-adically close to 1, so one might expect that less final precision is needed to bound the 
geometric Picard number. This expectation is confirmed by the following set of examples. 

Example 4.3.3. Table [2] lists some further instances in which the geometric Picard number 
condition in Proposition 14. 3. 11 can be verified (for p = 7, 11, 13, 17, 19), yielding the remaining 
cases of |^ Proposition 5.1]. We also list two cases where the geometric Picard number 
condition appears to hold but has not been verified definitively (for p = 23,29). In all 
cases, we took a = 1, (3 = 3, and final precision p^; the initial precision required was 
always p^, and for p < 19, the upper bound in (11. 5. 2. 21) was found to be 2. The examples 
for p = 11, 13, 17, 19 were found by prescreening with initial precision p^; the examples for 
p = 23, 29 were obtained by prescreening, but we did not complete the calculation at initial 
precision p'^. 



Table 2: Further instances and presumed instances of Proposition 14.3.11 over Fp. 



p 


/ 


CPU sec 


MB 


7 


—2x^ + 2x'^y + 2x'^w + y'^ + Sy'^w + yzw — yw"^ + 2z^ 


697 


63 


11 


—5x^ — 2x'^y — bxy"^ — 2xz'^ + — yz^ + 2z^ — Aw'^ 


5320 


106 


13 


3x^ — Qx'^z + y^ — 6yw'^ + 2z^ 


14997 


158 


17 


—x^ + 8x^y — xyw + y^ — y'^w + 2z^ + 5z^w — bzw"^ 


61996 


306 


19 


6x^ + 3x^z + 7xyw — 7xz^ + 8xzw — 9xw^ + y^ 


116323 


459 




—y^w — byz"^ + byw"^ + 2z^ — Az'^w — 2zw'^ 






23 


— llx'^ — dx'^y — 2x'^z — dx'^w — 3xyz 


? 


? 




—lOxyw + y^ + llyw"^ + 2z^ — Aw^ 






29 


4x^ + Ax'^w — 5x?/^ — lAxyw + y^ 


? 


? 




+7y'^z - 3yz'^ + 3yw'^ + 2z^ - 







4.4 Examples: degrees 4, 5 over F2 

To conclude, we edge towards the realm of coding theory proper, by considering examples 
over F2. Frustratingly, the asymptotic advantage obtained in the p-adic algorithms by taking 
p = 2 (which occur because the degree of the truncated Frobenius lift is a linear function of 
p) is somewhat counterbalanced by the need for additional precision to fight the unfortunate 
propensity of small integer denominators to be divisible by large powers of 2. Nonetheless, 
one can still say something. 

Example 4.4.1. The smooth quartic surface over F2 defined by the polynomial 

+ x^z + x^y^ + x^yw + x^z^ + x^zw + xy^ + y'^ + y^w + yz^ + z^ + z^w^ + 

was found to have Picard number 1 by computing a Frobenius matrix modulo 2^; this 
required initial precision 2^^, 7182 CPU seconds, and 472 MB of memory. This example was 
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found by prescreening using initial precision 2^, which required a mere 88 CPU seconds and 
52 MB of memory. 

If we pass from quartics to quintics, then middle cohomology becomes 53-dimensional, 
with Hodge numbers 4, 45, 4. (Again, the space we compute in is only 52-dimensional 
because passing from X to U removes one cycle class.) This removes the parity obstruction 
to having geometric Picard number 1, and one would expect to be able to find examples of 
such. Moreover, since in this case one is not forced to include the eigenvalue —2, which is 
indistinguishable from +2 modulo 2^, one might even expect to be able to work with less 
precision. The following example fulfills these expectations. 

Example 4.4.2. The smooth quintic surface over F2 defined by the polynomial 

+ x^yz + x'^y^w + x'^yz^ + x'^z^ + xz^w"^ 
+y^ + y^zw + y'^zw^ + yzw^ + z^ + z^w^ + 

was found to have geometric Picard number 1 by computing a Frobenius matrix modulo 2^; 
this required initial precision 2^^, 22685 CPU seconds, and 179 MB of memory. Note that in 
this case screening for the geometric Picard number is a nontrivial calculation, because we 
must check up to n = 210, which entails working in some large extensions of Q2. 

The example was found by prescreening with initial precision 2^. If we had needed final 
precision 2^^, we would have required initial precision 2^^; we project that in our implemen- 
tation, such a calculation would require roughly 100000 CPU seconds, and would have to be 
done in batches (or parallelized) to avoid memory overrun. 
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